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Abstract. We review the current cosmological status of neutrinos, with particular emphasis on their effects on Big Bang 
Nucleosynthesis, Large Scale Structure of the universe and Cosmic Microwave Background Radiation measurements. 
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INTRODUCTION 

The Standard Cosmological Model predicts the existence 
of a neutrino background (CvB) filling the universe with 
densities of the order n v w 1 10 cm -3 per flavor. Neutrino 
properties are rather difficult to be probed experimen- 
tally, due to the weakness of neutrino interactions which, 
especially at low energies, makes hopeless at present any 
perspective of direct detection of the CvB. Nevertheless, 
neutrinos are one of the most abundant relics of the pri- 
mordial universe and played a key role in different stages 
of its evolution. Several cosmological observables are 
then sensitive to neutrinos, and can be used to put bounds 
on their properties. 

Given their extremely low interaction rate, the natural 
out-of-equilibrium driving force of the expansion of the 
Universe pushed neutrinos to decouple from the thermal 
bath very early, when the temperature was 0{\ MeV). 
This temperature is close to the electron mass m e , setting 
the scale of the electron/positron annihilation, and both 
are close to the ^(0.1 MeV) scale of the synthesis of 
the light nuclei via thermonuclear fusion. So, Big Bang 
Nucleosynthesis (BBN) is a privileged laboratory for the 
CvB studies. In particular, it is sensitive to the v (weak) 
interactions as well as to the shape of the v e — V e phase 
space distributions entering the inter-conversion 
rates. Apart from the energy density due to the extra (i.e. 
non electromagnetic) relativistic degrees of freedom, the 
BBN tests the dynamical properties of the neutrinos in a 
thermalized (almost) CP-symmetric medium. 

Other cosmological probes are Cosmic Microwave 
Background (CMB) anisotropies or the Large Scale 
Structure (LSS) of the universe, which are, however, sen- 
sitive only to the CvB gravitational interaction. The role 
of neutrinos as the dark matter (DM) particles has been 
widely discussed since the early 1970s. For values of 
neutrino masses much larger than the present cosmic 



of the universe, cosmology 



temperature one finds a contribution in terms of the crit- 
ical density H v ~ 0.0108 /i~ 2 £ra v /eV, h being the Hub- 
ble parameter in units of 100 Km s _1 Mpc -1 . Nowa- 
days, we know that neutrinos cannot constitute all the 
DM (£2dm ~ dm ~ 0-3 d), and the main question is 
how large the contribution of neutrinos can be, deduc- 
ing £2 V from their contribution to cosmological pertur- 
bations. In fact, neutrino background erases the density 
contrasts on wavelengths smaller than a mass-dependent 
free-streaming scale. Neutrinos of sub-eV mass behave 
almost like a relativistic species for CMB considerations 
and therefore the power spectrum suppression can be 
seen only in LSS data. Even if neutrino mass influences 
only slightly the spectrum of CMB anisotropies, it is cru- 
cial to combine CMB and LSS observations, because the 
former give independent constraints on the cosmological 
parameters, and partially remove the parameter degener- 
acy that would arise in an analysis of the LSS only. 

NEUTRINOS AND COSMOLOGY 

At temperatures above 0(\ MeV), neutrinos are in ther- 
mal equilibrium with the thermal bath and their distribu- 
tion is a perfect Fermi-Dirac one, 

where y = p/T v and E, a = fl a /T v (here jj. a is the chem- 
ical potential of the flavor 0£, which is neglected in 
the standard scenario). As the temperature goes down, 
the universe expansion prevents weak interactions from 
maintaining neutrinos in equilibrium and they decou- 
ple. As a first approximation, the neutrino decoupling 
can be described as an instantaneous process taking 
place around 2-4 MeV, without any overlap in time with 
e + — e~ annihilation. All flavors would then keep Fermi- 



Dirac distributions (both neutrino momenta and temper- 
ature redshift identically with the universe expansion), 
but the neutrino temperature T v will not benefit of the 
entropy release from e + — e~ annihilations. The asymp- 
totic ratio T / T v for T <C m e can be evaluated in an ana- 
lytic way, and turns out to be ( 1 1 /4) 1 / 3 ss 1 .40 1 . 

More accurate calculations by solving the kinetic 
equations have been performed, and they show a par- 
tial entropy transfer to the neutrino plasma. As a conse- 
quence, the neutrino distributions get distorted, since this 
transfer is more efficient for larger neutrino momenta. In 
HI it was shown that with a very good approximation 
the distortion in the a-th flavor can be described as 
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with x = m e /T v . The electron neutrinos get a larger en- 
tropy transfer than the jx and T, since they also interact 
via charged currents with the e plasma. The effective 
ratio T /T v « 1.3984 is slightly lower than the instanta- 
neous decoupling estimate. 

The incomplete decoupling of neutrinos also induces 
a modification in the contribution of neutrinos to the 
energy density. By fully consistently including order a 
QED corrections to the photon and e ± equation of state, 
in 01 the energy density in the neutrino fluid is found 
to be enhanced by 0.935% (for v e ) and 0.390% (for v M 
and v T ). A refined treatment, also including the effects 
of three-flavor neutrino oscillations, has been recently 
provided in |4]. 

The contribution of neutrinos to the total relativistic 
energy density of the universe is usually parameterized 
via the "effective number" of neutrinos, N e ff, 



N, 



eff 



Py. 



(3) 



N e ff measures neutrino energy density in "units" of the 
energy density of massless neutrinos with zero chemi- 
cal potential, but it can in principle receive a contribu- 
tion from other (relativistic) relics. For three massless 
neutrinos with zero chemical potential and in the limit 
of instantaneous decoupling, N e ff = 3. The inclusion of 
entropy transfer between neutrinos and the thermal bath 
modifies this number to about 3.04 at the CMB epoch. 

Shortly after neutrino decoupling the temperature 
reaches the value of the neutron-proton mass difference, 
and weak interactions are no longer fast enough to main- 
tain equilibrium among nucleons: a substantial final neu- 
tron fraction survives, however, down to the phase of 
nucleosynthesis where all neutrons become practically 
bound in 4 He nuclei. The predicted value of the 4 He mass 
abundance, Y p , is poorly sensitive to the nuclear network 
details and has only a weak, logarithmic dependence on 
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FIGURE 1. The relative correction to the n — » p (solid line) 
and p — > n (dashed line) total rates, due to the neutrino distor- 
tion (see Ref. 1 7] for details). 



the baryon fraction of the universe, COj, = £2f,/z 2 , being 
fixed essentially by the ratio of neutron to proton number 
density at the onset of nucleosynthesis. This in turn cru- 
cially depends on the weak rates and on the (standard or 
exotic) neutrino properties. 

In several papers (see 0101 and references therein) the 
value of Y p has been computed by improving the evalua- 
tion of the weak rates including electromagnetic radiative 
corrections, finite mass corrections and weak magnetism 
effects, as well as the plasma and thermal radiative ef- 
fects. In particular in |7] it has been also considered the 
effect of the neutrino spectra distortions and of the pro- 
cess 7+ p <-> v e + e + + n, which is kinematically forbid- 
den in vacuum, but allowed in the thermal bath. The lat- 
ter is shown to give a negligible contribution, while neu- 
trino distortions have a significant influence on the rates 
for different reasons: a) the larger mean energy of v e in- 
duces a 8T Ve and (indirectly, through a decrease in p e m ) 
a 8T; b) the ratio T/T v , which enters the weak rates, is 
changed. Moreover, the time-temperature relationship is 
changed and so the time at which BBN starts. 

The total effect on the rates are shown in Figure ^ 
Even though one would expect effects up to 0(\%), 
the spectral distortion and the changes in the energy 
density and T V (T) conspire to almost cancel each other, 
so that Y p is changed by a sub-leading 0(0.1%). This 
effect is of the same order of the predicted uncertainty 
coming from the error on the measured neutron lifetime, 
T„ = 885.7 ± 0.8 s 1 8], and has to be included in quoting 
the theoretical prediction, Y p = 0.248 1 ± 0.0004 ( 1 a, for 
co b = 0.023 ±0.001). 

Apart from the uncertainty on a>b, the 2 H, 3 He and 7 Li 
abundance predictions are mainly affected by the nuclear 



reaction uncertainties. An updated and critical review of 
the nuclear network and a new protocol to perform the 
nuclear data regression has been widely discussed in 0], 
to which we address for details. 

A different scenario, the Degenerate BBN (DBNN) 
one, has received a new attention in the last few years, 
especially when the first data on CMB seemed to indicate 
a tension between the determination of ©/, from CMB 
and standard BBN |9]. Such a tension could, in fact, 
be relaxed assuming that the number of v a and v a be 
different, i.e. jx a ^ in Eq.([l). In this degenerate case, 
changes are expected: a) in the weak rates (a reduction 
in Y p ), since a positive £, e enhances n — > p processes 
with respect to the inverse ones, and modifies the initial 
condition on the n/p ratio at J > 1 MeV; b) in the 
expansion rate (an increase in Y p ), since non zero ^ a 's 
contribute to total N e g- as 
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When several neutrino species are degenerate, both the 
effects might combine and particular values of ^ a 's ex- 
ist for which the predictions of DBBN are still in good 
agreement with the observational data on the abundances 
of primordial elements. Notice that BBN is more sensi- 
tive to neutrino degeneracy than CMB or LSS, due to the 
further effect played by the v e — V e distributions in the 
weak rates, while the latter are only sensitive to the extra 
energy density present in the | ^ case. 

Earlier claims of discrepancies between BBN and 
CMB have been largely overcome by new data, and it 
was recently realized 1 10] that the flavor oscillations in 
the primordial plasma induced by (presently determined) 
mass differences and mixing angles from atmospheric 
and solar neutrinos almost equalize the three asymme- 
try parameters Still, exotic models, where both a 
common relatively large £, and a N e ff 7^ 3.04 exist, have 
been considered ("hidden relativistic degrees of free- 
dom" 111 ll ll2IO . and were shown to be compatible with 
the data. However, if one sticks to the scenario with N e ff 
fixed by standard Physics, introduces no sterile species, 
and assumes the CMB prior on ©/,, BBN turns into a 
poweful "leptometer", constraining the common t, to an 
unprecedent accuracy even under conservative assump- 
tions for Y p ltl3ll . This provides an indirect consistency 
check for the sphaleron mechanism at electroweak phase 
transition, predicting baryon and lepton asymmetries of 
the same order. 

Let us come to neutrino role in structure formation. 
In general, neutrinos tend to stream freely across grav- 
itational potential wells, and to erase density perturba- 
tions. Free-streaming is efficient on a characteristic scale 
Aj called the Jeans length, corresponding roughly to the 
distance on which neutrinos can travel in a Hubble time. 



For ultra-relativistic neutrinos, Aj is by definition equal 
to the Hubble radius c/H, but for non-relativistic ones it 
is lower than c/H. Neutrinos with masses smaller than 
approximately 0.3 eV are still relativistic at the time of 
last scattering, and their direct effect on the CMB per- 
turbations is identical to that of massless neutrinos. In 
the intermediate mass range from 10~ 3 eV to 0.3 eV, the 
transition to the non-relativistic regime takes place dur- 
ing structure formation, and the matter power spectrum 
will be directly affected in a mass-dependent way. Wave- 
lengths A smaller than the current value of the neutrino 
Aj are suppressed by free-streaming. The largest wave- 
lengths, which remain always larger than the neutrino Aj, 
are not affected. Finally, there is a range of intermediate 
A which become smaller than the neutrino Aj for some 
time, and then encompass it again: these scales smoothly 
interpolate between the two regimes. The net signature 
in the matter power spectrum P(X , t) is a damping of all 
wavelengths smaller than the Hubble scale at the time To 
of the transition of neutrinos to a non-relativistic regime. 

Then, for A <C Aj, if Q. v <C £2 m the suppression is given 
roughly by the factor 
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that is by the ratio between neutrino and matter energy 
densities. 

Notice that one can somehow play with both N e g and 
Y* m v and find models which give excellent fits of the 
data. In fact, models with massive neutrinos have sup- 
pressed power at small scale. Adding relativistic energy 
further suppresses power at scales smaller than the hori- 
zon at matter-radiation equality. For the same matter den- 
sity such a model would therefore be even more incom- 
patible with data. However, if the matter density is in- 
creased together with m v and N e ff, data can be described 
very nicely (see, for example Figure 3 in 1 14]). 



COMPARISON WITH DATA 

Still few years ago, the BBN theory together with the 
observations of the abundances of primordial nuclides 
were used to determine the baryon fraction of the uni- 
verse, (Ob- Nowadays (Oy ~ 0.023 is fixed to better than 
5% accuracy by detailed CMB anisotropics analysis 1 1], 
thus leaving the BBN as an over-constrained and very 
predictive theory. Once cot, = 0.023 ±0.001 is plugged 
into the BBN theory, the prediction for the deuterium, 
which is the nuclide most sensitive to C0b, nicely fits the 
range of the observed values in high redshift, damped 
Ly-a QSO systems 1 15], thus offering a remarkable ex- 
ample of internal consistency of the current cosmological 
scenario. Moreover, the predictions of other light nuclei 



TABLE 1. Bounds on N e ff (2 o") from different analyses 



TABLE 3. Bounds on Y* m v (2 o") from different analyses 



Ref. 


Bound on N e ff 


Data used 


Ljj 


l-8<N eff <3.7 


CMB,BBN 


[12] 


1.3 <N eff <6.1 


CMB, BBN(D) 


[12] 


1.6<N eff <3.6 


BBN(D+7 p ) 


[17] 


1.4<N eff <6.8 


CMB, LSS, HST 


[18] 


1.9(2.3) < N eff < 7.0(3.0) 


CMB, LSS, (+BBN) 


[19] 


1.7<N eff <3.0 


CMB, BBN 


[20] 


N eff < 4.6 


CMB, BBN 


[21] 


1.90 < N eff < 6.62 


CMB, LSS, HST 



Ref. 


Bound on £m v 


Data used 


UJ 


<0.69 


CMB, LSS, Lya 


[18] 


< 1.01 


CMB, LSS, HST, SNIa 


[22] 


<0.65 


CMB, LSS, HST, Lya 


[23] 


- 56+ 030 


CMB, LSS, i gas , XLF 


[24] 


< 1.7 


CMB, LSS 


[25] 


<0.75 


CMB, LSS, HST 


[26] 


< 1.0(0.6) 


CMB, LSS (+HST, SNIa) 


[27] 


<0.47 


CMB, LSS, HST, SNIa, Lya 



which at least qualitatively agree with the observed val- 
ues are likely to put constraints on the Galactic chemical 
evolution ( 3 He) or on the temperature scale calibration or 
depletion mechanisms in PopII halo stars ( Li). 

The determination of Y p is usually performed by ex- 
trapolating to zero metallicity the measurements done in 
dwarf irregular and blue compact galaxies. The typical 
statistical errors are of the order of 0.002 {i.e., at the 1% 
level), but the systematics are such that in the recent re- 
analysis flrjll the authors argue for the conservative range 
0.232 <Y p < 0.258, i.e. a 1 a error of ff(5%). 

In Table ^ tne present bounds on the effective num- 
ber of neutrinos from various analyses are presented, to- 
gether with the type of data used. The most stringent 
bounds come from BBN alone (Deuterium+Helium), 
while CMB and BBN-Deuterium are less effective in 
constraining N e ff. Some differences are due to (slightly) 
different databases or assumptions. 

Table|2shows the bounds on <jj e . The interval from Ref. 
lll2ll is broader than what in Ref. full , since in the first 
case only a prior from CMB instead of all data is used in 
the analysis. The third line shows the bound obtained in 
lll3ll assuming only standard physics, while the previous 
two bounds assume no prior on N e ff. 

Finally, Table [5] shows, with the same notations, the 
upper bound on the neutrino mass. As can be gauged 
from this table, a fairly robust bound on the sum of 
neutrino masses is at present somewhere around 1 eV, 
depending on the specific priors and data sets used. 

In conclusion, the complementarity among different 
fields of cosmology (BBN, CMB, LSS) can be used to 
test the role of neutrinos from very early epochs (redshift 
Z ~ 10 10 ) down to relatively recent history (z ~ a few). 

TABLE 2. Bounds on t, e (2 a) from different analyses 
Ref. Bound on % e Data used 

|T1J -0. 10 < t, e < 0.25 CMB, BBN 

[121 -0.13 < ke< 0.31 BBN(D+y p )+prior on co b 

[13] -0.05 < 4. < 0.07 BBN(y^)+priors on (O h , N eff 



New nuclear rate measurements and a better under- 
standing of possible systematics affecting primordial 
abundance determination (BBN) and more data together 
with an increased precision (CMB and LSS) will give us 
the opportunity to constrain standard neutrino properties 
as well as to test new physics in the neutrino sector, 
gaining at the same time a deeper insight on the physics 
of the early universe. 
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